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Single  photon  detection  is  a  requisite  technique  in  quantum-optics  experiments 
in  both  the  optical  and  the  microwave  domains.  However,  the  energy  of  mi¬ 
crowave  quanta  are  four  to  five  orders  of  magnitude  less  than  their  optical 
counterpart,  making  the  efficient  detection  of  single  microwave  photons  ex¬ 
tremely  challenging.  Here,  we  demonstrate  the  detection  of  a  single  microwave 
photon  propagating  through  a  waveguide.  The  detector  is  implemented  with 
an  “impedance-matched”  artificial  A-system  comprising  the  dressed  states  of 
a  driven  superconducting  qubit  coupled  to  a  microwave  resonator.  We  attain  a 
single-photon  detection  efficiency  of  0.66  ±  0.06  with  a  reset  time  of  ~  400  ns. 
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This  detector  can  be  exploited  for  various  applications  in  quantum  sensing, 
quantum  communication  and  quantum  information  processing. 

Single -photon  detection  is  essential  to  many  quantum-optics  experiments,  enabling  photon 
counting  and  its  statistical  and  correlational  analyses  (/).  It  is  also  an  indispensable  tool  in 
many  protocols  for  quantum  communication  and  quantum  information  processing  (2-5).  In  the 
optical  domain,  various  kinds  of  single-photon  detectors  are  commercially  available  and  com¬ 
monly  used  ( 1 ,  6).  However,  despite  the  latest  developments  in  nearly-quantum- limited  ampli¬ 
fication  (7,  8)  and  homodyne  measurement  for  extracting  microwave  photon  statistics  (9),  the 
detection  of  a  single  microwave  photon  in  an  itinerant  mode  remains  a  challenging  task  due  to 
its  correspondingly  small  energy.  Meanwhile,  the  demand  for  such  detectors  is  rapidly  increas¬ 
ing,  driven  by  applications  involving  both  microwave  and  hybrid  optical-microwave  quantum 
systems. 

In  this  report  we  demonstrate  an  efficient  and  practical  single-microwave-photon  detector 
based  on  the  deterministic  switching  in  an  artificial  A-type  three-level  system  implemented 
using  the  dressed  states  of  a  driven  superconducting  quantum  circuit.  The  detector  operates 
in  a  time-gated  mode  and  features  a  high  quantum  efficiency  0.66  ±  0.06,  a  low  dark-count 
probability  0.014  ±  0.001,  a  bandwidth  ~  2tt  x  8  MHz,  and  a  fast  reset  time  ~  400  ns.  It  can 
be  readily  integrated  with  other  components  for  microwave  quantum  optics. 

Our  detection  scheme  carries  several  advantages  compared  with  previous  proposals.  It  uses 
coherent  quantum  dynamics,  which  minimizes  energy  dissipation  upon  detection  and  allows 
for  rapid  resetting  with  a  resonant  drive,  in  contrast  to  schemes  that  involve  switching  from 
metastable  states  of  a  current-biased  Josephson  junction  into  the  finite  voltage  state  (10-12). 
Moreover,  our  detection  scheme  does  not  require  any  temporal  shaping  of  the  input  photons, 
nor  precise  time-dependent  control  of  system  parameters  adapted  to  the  temporal  mode  of  the 
input  photons,  in  contrast  to  recent  photon-capturing  experiments  (13-15).  It  also  achieves  a 
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high  efficiency  without  cascading  many  devices  in  series  (10, 16). 

The  operating  principle  of  the  detector  fully  employs  the  elegance  of  waveguide  quantum 
electrodynamics,  which  has  recently  attracted  significant  attention  in  various  contexts  surround¬ 
ing  photonic  quantum  information  processing  (17-20).  When  electromagnetic  waves  are  con¬ 
fined  and  propagate  in  an  one-dimensional  (ID)  mode,  their  interaction  with  a  quantum  emit- 
ter/scatterer  is  substantially  simplified  and  enhanced  compared  with  three-dimensional  modes. 
These  advantages  result  from  the  natural  spatial-mode  matching  of  the  emitter/scatterer  with 
a  ID  mode  and  its  resulting  enhancement  of  quantum  interference  effects.  Remarkable  exam¬ 
ples  are  the  perfect  extinction  of  microwave  transmission  for  an  artificial  atom  coupled  to  a  ID 
transmission  line  (21,  22),  the  photon-mediated  interaction  between  two  remote  atoms  coupled 
to  a  ID  transmission  line  (23),  and  the  perfect  absorption  -  and  thus  “impedance  matching”  -  of 
a  A-type  three-level  system  terminating  a  ID  transmission  line  (24, 25).  In  the  latter  system,  the 
incident  photon  deterministically  induces  a  Raman  transition  which  changes  the  state  of  the  A 
system  (24, 26).  This  effect  has  recently  been  demonstrated  in  both  the  microwave  and  optical 
domains  (25,  27),  indicating  its  potential  for  photon  detection  (28)  as  well  as  for  implementing 
deterministic  entangling  gates  with  photonic  qubits  (29). 

Our  device  consists  of  a  superconducting  flux  qubit  capacitively  and  dispersively  coupled 
to  a  microwave  resonator  (Fig.  IB)  (30).  With  a  proper  choice  of  the  qubit  drive  frequency  oj,\ 
and  power  Pd,  the  system  functions  as  an  impedance-matched  A-system  with  identical  radiative 
decay  rates  from  its  upper  state  to  its  two  lower  states  (Fig.  1A)  (24,  25).  The  qubit-resonator 
coupled  system  is  connected  to  a  parametric  phase-locked  oscillator  (PPLO),  which  enables  fast 
and  non-destructive  qubit  readout  (31). 

Figure  1C  shows  the  level  structure  of  the  qubit-resonator  system  and  the  protocol  for  the 
single  photon  detection.  We  label  the  energy  levels  |  q,n)  and  their  eigenfrequencies  u\q,n), 
where  q  =  {g,  e}  and  n  —  {0, 1,  •  •  •  }  respectively  denote  the  qubit  state  and  the  photon  number 
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in  the  resonator.  In  the  dispersive  coupling  regime,  the  qubit-resonator  interaction  renormalizes 
the  eigenfrequencies  to  yield  uj\g,n)  =  nuJr  and  u\e,n)  =  ^'ge  +  n( uT  —  2%),  where  c uge  and  ur  are 
the  renormalized  frequencies  of  the  qubit  and  the  resonator,  respectively,  and  x  is  the  dispersive 
frequency  shift  of  the  resonator  due  to  its  interaction  with  the  qubit.  Only  the  lowest  four  levels 
with  n  =  0  or  1  are  relevant  here. 

We  prepare  the  system  in  its  ground  state  | g,  0)  (Fig.  1C,  Initialization)  and  apply  a  drive 
pulse  to  the  qubit  (Fig.  1C,  Detection).  In  a  frame  rotating  at  cud,  the  level  structure  becomes 
nested,  i.e.,  cu|Si0)  <  w|e,o>  <  ^le.i)  <  u\g,i)^  f°r  in  the  range  coge  —  2\  <  cud  <  cuge.  On  the 
plateau  of  the  drive  pulse  (Pd  >  0),  the  lower-two  levels  \g,  0)  and  |e,  0)  (high-two  levels  | g,  1) 
and  |e,  1))  hybridize  to  form  dressed  states  |1)  and  |2)  (|3)  and  |4)).  Under  a  proper  choice  of  Pd, 
the  two  radiative  decay  rates  from  |4)  (or  |3))  to  the  lowest  two  levels  become  identical.  Thus, 
an  impedance-matched  A  system  comprising  1 1),  |2),  and  |4)  (alternatively,  |1),  |2),  and  |3))  is 
realized.  An  incident  single  microwave  photon  (Gaussian  envelope,  length  ts ),  applied  through 
the  signal  port  and  in  resonance  with  the  1 1)  — >  |4)  transition,  deterministically  induces  a  Raman 
transition,  |1)  — >  |4)  — >•  |2),  and  is  down-converted  to  a  photon  at  the  |4)  — >  |2)  transition 
frequency.  This  process  is  necessarily  accompanied  by  an  excitation  of  the  qubit  (24, 25). 

To  detect  the  photon,  we  adiabatically  switch  off  the  qubit  drive  and  dispersively  read  out 
the  qubit  state  (Fig.  1C,  Readout).  We  apply  a  readout  pulse  with  the  frequency  turd  =  uT  —  2%  = 
cu|eii)  —  cu|e)0)  through  the  signal  port,  which,  upon  reflection  at  the  resonator,  acquires  a  qubit- 
state-dependent  phase  shift  of  0  or  n.  This  phase  shift  is  detected  by  the  PPLO  with  high 
fidelity  (31).  In  the  present  setup,  the  readout  fidelity  of  the  qubit  is  ~  0.9,  which  is  primarily 
limited  by  qubit  relaxation  prior  to  readout  (31). 

We  first  determine  the  operating  point  where  the  A  system  deterministically  absorbs  a  signal 
photon.  We  simultaneously  apply  a  drive  pulse  of  length  td  =  178  ns  and  a  signal  pulse  of 
length  ts  =  85  ns,  and  proceed  to  measure  the  reflection  coefficient  |r|  of  the  signal  pulse  as 
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a  function  of  the  drive  power  Pd  and  the  signal  frequency  uB  (Fig.  2A).  The  signal  pulse  is  a 
weak  coherent  state  with  mean  photon  number  ns  ~  0.1.  A  pronounced  dip  with  a  depth  of 
<—25  dB  is  observed  in  |r|  at  ( Pd,Lus/2n )  =  (—76  dBm,  10.268  GHz),  in  close  agreement 
with  theory  (Fig.  2C).  The  dip  indicates  a  near-perfect  absorption  condition,  i.e.  impedance 
matching,  where  the  reflection  of  the  input  microwave  photon  vanishes  due  to  destructive  self¬ 
interference.  Correspondingly,  a  deterministic  Raman  transition  of  1 1)  — >■  |4)  — »  |2)  is  induced, 
and  the  qubit  state  is  flipped. 

To  obtain  a  ‘click’  corresponding  to  single-photon  detection,  we  read  out  the  qubit  state  by 
using  the  PPLO  immediately  after  the  Raman  transition.  Before  initiating  readout,  the  drive 
pulse  is  turned  off  to  suppress  unwanted  Raman  transitions  induced  by  the  readout  pulse,  e.g., 
|2)  — »  |3)  — >  |1).  We  repeatedly  apply  the  pulse  sequence  in  Fig.  1C  104  times  and  evaluate 
the  single-photon  detection  efficiency  q  =  P(|e))/[1  —  P(0)],  where  P(  |e))  and  P(0)  are  the 
probabilities  for  the  qubit  being  in  the  excited  state  and  the  signal  pulse  being  in  the  vacuum 
state,  respectively.  Figure  2B  depicts  q  as  a  function  of  Pd  and  us.  The  dark  count  probability 
of  the  detector  —  mainly  caused  by  the  nonadiabatic  qubit  excitation  due  to  the  drive  pulse 
and  the  imperfect  initialization  —  is  subtracted  when  evaluating  q  (32).  We  observe  that  q  is 
maximized  at  the  dip  position  in  Fig.  2 A  in  accordance  with  the  impedance-matching  condition. 
We  also  confirm  that  the  result  agrees  with  numerical  calculations  based  on  the  parameters 
determined  independently  (Fig.  2D).  The  maximum  value,  q  =  0.66  ±  0.06,  is  obtained  at 
(Pd,tus/27r)  =  (—75.5  dBm,  10.268  GHz)  (Fig.  2E)  (32).  The  efficiency  exceeds  0.5  over  a 
signal-frequency  range  of  ~  20  MHz,  which  is  comparable  to  the  bandwidth  of  the  detector, 
k/2i r  ~  16  MHz  (32). 

In  the  Fig.  3A,  we  plot  efficiency  q  as  a  function  of  the  signal  pulse  length  ts.  Here,  we 
fix  cos  and  Pd  at  the  values  which  maximize  q  in  Fig.  2E.  The  drive  pulse  duration  t,\  is  set 
to  be  fd  =  1.5ts  +  50  ns,  which  empirically  maximizes  q  at  each  ts.  We  observe  that  q  is  a 
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non-monotonic  function  of  ts  and  attains  a  maximum  at  ts  =  85  ns.  The  initial  increase  of  77  at 
short  ts  is  due  to  the  narrowing  of  the  signal  bandwidth  resulting  in  an  improved  overlap  with 
the  detection  bandwidth.  For  longer  ts,  the  qubit  relaxation  limits  77  (28).  Next,  we  examine 
how  the  photon  detector  behaves  when  ns  in  the  signal  pulse  is  varied.  Figure  3B  shows  77  as  a 
function  of  ns  for  fixed  signal  pulse  lengths  at  ts  =  34,  85,  and  189  ns.  The  detection  efficiencies 
stay  constant  for  ns  <  1  regardless  of  the  pulse  lengths.  This  validates  the  determination  of  77 
in  our  measurements  using  signal  pulses  in  the  weak  coherent  states.  For  ns  >  1,  77  slightly 
depends  on  ns  because  of  the  possibility  to  drive  multiple  Raman  transitions. 

After  a  single-photon  detection  event,  the  qubit  remains  in  the  excited  state  until  it  sponta¬ 
neously  relaxes  to  the  ground  state,  which  leads  to  a  relatively  long  dead  time  of  the  detector. 
However,  our  coherent  approach  allows  us  to  implement  a  fast  reset  protocol  (Fig.  4A):  in  con¬ 
junction  with  the  drive  pulse  that  forms  the  A  system,  we  apply  a  relatively  strong  reset  pulse 
(signal  port)  which  induces  an  inverse  Raman  transition,  |2)  — >  |3)  — >  |1).  We  optimize  the 
drive-pulse  power  Pdr  and  the  reset-pulse  frequency  u>rst  such  that  the  resulting  qubit  excita¬ 
tion  probability  P(|e))  is  minimized  (Fig.  4B).  At  the  optimal  reset  point  (Pi,., u7rsti/27r)  = 
(—72.1  dBm,  10.162  GHz),  P(|e))  attains  a  minimum  value  0.017  ±  0.002,  equivalent  to  the 
value  0.016  ±  0.001  obtained  in  the  absence  of  the  initial  7r -pulse  used  to  mimic  a  photon  ab¬ 
sorption  event.  Without  a  reset  pulse,  we  obtain  P(|e))  =  0.490  ±  0.010.  A  comparison  of  the 
two  results  indicates  that  the  reset  pulse  is  highly  efficient.  Moreover,  we  confirm  that  the  reset 
protocol  does  not  affect  the  succeeding  detection  efficiency  and  that  the  time-gated  operation 
can  be  repeated  in  a  rate  exceeding  1  MHz  (32). 

For  the  moment,  the  detection  efficiency  of  this  detector  is  limited  by  the  relatively  short 
qubit  relaxation  time,  T\  ~  0.7  /is.  Nonetheless,  our  theoretical  work  indicates  that  efficien¬ 
cies  reaching  ~  0.9  are  readily  achievable  with  only  a  modest  improvement  of  the  qubit  life¬ 
time  (28).  An  extension  from  the  time-gated-mode  to  the  continuous-mode  operation  is  also 
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possible  (55). 
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Figure  1 :  Experimental  setup  and  pulse  sequence.  (A)  Image  of  the  sample  chip  containing 
a  flux  qubit  and  a  superconducting  microwave  resonator  coupled  capacitively  and  operated  in 
the  dispersive  regime.  For  certain  proper  conditions  of  the  qubit  drive,  the  coupled  system 
functions  as  an  impedance-matched  A-type  three-level  system.  (B)  Schematic  of  the  itinerant 
microwave-photon  detector  consisting  of  the  coupled  system  and  connected  to  a  parametric 
phase-locked  oscillator  (PPLO)  via  three  circulators  in  series.  The  circuit  has  three  input  ports: 
signal,  qubit  drive,  and  pump  for  the  PPLO.  (C)  Energy-level  diagram  of  the  coupled  system 
and  the  pulse  sequence  for  single-photon  detection.  The  system  is  first  prepared  in  the  ground 
state.  During  the  detection  stage,  we  concurrently  apply  the  drive  and  signal  pulses.  The  drive 
is  parameterized  to  fulfill  the  impedance-matched  condition  such  that  a  signal  photon  (blue 
arrow)  induces  a  deterministic  Raman  transition.  A  down-converted  photon  (green  arrow)  is 
emitted  in  the  process  and  discarded.  In  the  readout  stage,  we  detect  the  qubit  excited  state 
nondestructively  by  sending  a  qubit  readout  pulse.  The  qubit-state-dependent  phase  shift  in  the 
reflected  pulse  is  discriminated  by  the  PPLO. 
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Figure  2:  Impedance  matching  and  itinerant  microwave-photon  detection.  (A)  Amplitude 
of  the  reflection  coefficient  |r|  of  the  input  signal  pulse  with  mean  photon  number  ns  ~  0.1  as 
a  function  of  the  qubit  drive  power  Pd  and  the  signal  frequency  cos.  The  PPLO  is  not  activated 
during  this  measurement.  The  impedance-matched  point  is  resolved  (dark-blue  region),  where 
the  input  microwave  photon  is  absorbed  almost  completely.  In  the  inset,  we  also  observe  another 
dip  in  \r\,  corresponding  to  the  Raman  transition  of  1 1)  — >■  |3)  — *  |2).  Microwave  power  levels 
stated  in  this  article  are  referred  to  the  value  at  the  corresponding  ports  on  the  sample  chip. 

(B)  Detection  efficiency  //  of  an  itinerant  microwave  photon.  The  efficiency  hits  the  maximum 
at  the  impedance-matched  point,  where  the  Raman  transition  of  1 1)  — >  |4)  — >•  \2)  takes  place. 

(C)  and  (D)  Theoretical  predictions  corresponding  to  A  and  B.  (E)  Cross-sections  of  B  (blue 
dots)  and  D  (red  dashed  line)  at  uis/2tt  =  10.268  GHz.  The  error  bars  are  due  to  the  uncertainty 
in  the  input  power  calibration. 
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Figure  3:  Optimization  of  the  efficiency.  (A)  Single  photon  detection  efficiency  rj  as  a  function 
of  the  signal  pulse  length  ts.  The  mean  photon  number  ns  for  the  weak-coherent  signal  pulse  is 
~  0.1.  (B)  q  as  a  function  of  ns.  Dashed  lines  indicate  theoretical  predictions. 
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Figure  4:  Demonstration  of  the  fast  reset  protocol.  (A)  Pulse  sequence  used  to  evaluate  the 
reset  efficiency.  The  initial  7r -pulse  mimics  a  single-photon  detection  and  excites  the  qubit. 
During  the  reset  stage,  a  drive  pulse  and  a  reset  pulse  with  the  mean  photon  number  of  nrst  ~ 
43  are  concurrently  applied,  inducing  an  inverse  Raman  transition:  |2)  — >  |3)  — >  |1).  The 
remaining  population  in  the  |e)  state  is  then  detected.  (B)  Population  of  the  qubit  excited  state 
after  the  reset  operation  P(|e)),  as  a  function  of  the  reset-pulse  frequency  ujIst  and  the  drive- 
pulse  power  Pdr.  (C)  Theoretical  prediction  for  B  with  no  free  parameters.  (D)  Cross  sections 
of  B  (blue  dots)  and  C  (red  dashed  line)  at  curst/27r  =  10.162  GHz. 
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Device 


Our  device  is  composed  of  a  A/2  superconducting  coplanar  waveguide  (CPW)  resonator  and 
a  superconducting  flux  qubit  (Fig.  1A).  The  CPW  resonator  is  made  of  a  50-nm-thick  Nb  film 
sputtered  on  a  300-//m-thick  undoped  silicon  wafer  with  a  300-nm-thick  thermal  oxide  on  the 
surface.  It  is  patterned  by  electron-beam  (EB)  lithography  using  the  ZEP520A-7  resist  and 
CF4  reactive  ion  etching.  The  flux  qubit  with  three  Josephson  junctions,  where  one  is  made 
smaller  than  the  other  two  by  a  factor  of  a,  is  fabricated  by  EB  lithography  and  double-angle 
evaporation  of  A1  using  PMMA  (50  nm)/Ge  (50  nm)/MMA  (400  nm)  trilayer  resist  (Fig.  SIB). 
The  thicknesses  of  the  bottom  and  the  top  A1  layers  separated  by  an  AI2O3  layer  are  20  and  30 
nm,  respectively.  In  order  to  realize  a  superconducting  contact  between  Nb  and  Al,  the  surface 
of  Nb  is  cleaned  by  Ar  ion  milling  before  the  evaporation  of  Al.  The  qubit  is  located  at  one  end 
of  the  resonator  and  is  coupled  to  the  resonator  dispersively  through  a  capacitance  of  4  fF,  while 
it  is  coupled  to  the  drive  port  inductively  (Fig.  SI  A). 

The  flux  qubit  is  always  biased  with  a  half  flux  quantum  where  the  transition  frequency  of 
the  qubit  cuge  from  the  ground  state  | g)  to  the  excited  state  \e)  is  2tt  x  5.508  GHz  (7\  ~  700  ns 
during  photon-detection  experiments),  while  the  resonator  frequency  cer  is  27t  x  10.256  GHz 
( Q  factor  ~  630)  when  the  qubit  is  in  the  | g)  state.  It  is  shifted  by  a  dispersive  interaction  with 
the  qubit  of  —  2\  =  —2n  x  69  MHz,  which  is  enhanced  by  the  straddling  effect  (i)  and  the 
capacitive  coupling  (2)  when  the  qubit  is  in  the  \e)  state.  Note  that  uge  and  ur  denote  not  their 
bare  frequencies  but  the  renormalized  ones  including  the  dispersive  shifts  (3). 

A  parametric  phase-locked  oscillator  (PPLO)  (4),  which  is  previously  operated  as  a  flux- 
driven  Josephson  parametric  amplifier  (JPA)  (5)  consists  of  a  A/4  superconducting  CPW  res¬ 
onator  terminated  by  a  dc-SQUID  (superconducting  quantum  interference  device).  A  pump 
port  is  coupled  to  the  SQUID  loop  inductively.  The  device  was  fabricated  by  the  planarized  Nb 


Figure  SI:  Image  of  the  qubit-resonator  coupled  system.  (A)  False-colored  optical  image  of 
the  device  magnified  at  the  qubit  part.  The  qubit  (white)  is  coupled  to  the  center  conductor  of 
the  coplanar  waveguide  (CPW)  resonator  (green)  through  a  capacitance  of  4  fF.  (B)  Scanning 
electron  micrograph  of  the  three -junction  flux  qubit.  The  areas  shaded  by  yellow  indicate  the 
Josephson  junctions. 
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trilayer  process  at  MIT  Lincoln  Laboratory.  The  resonator  and  the  pump  port  are  made  out  of 
a  150-nm-thick  Nb  film  sputtered  on  a  Si  substrate  covered  by  a  500-nm-thick  Si02  layer.  A 
static  resonant  frequency  of  the  PPLO  is  c cPO  =  27 r  x  10.948  GHz.  The  PPLO  chip  is  the  same 
as  the  one  used  in  Ref.  (4). 

Experimental  setup 

A  schematic  of  the  measurement  setup  including  the  wiring  in  a  cryogen-free  3He/4He  dilution 
refrigerator,  circuit  components,  and  instruments  used  in  the  experiment  is  shown  in  Fig.  S2. 

The  qubit+resonator  and  the  PPLO  circuits  are  fabricated  on  separate  chips  and  are  sep¬ 
arately  mounted  in  microwave-tight  packages  equipped  with  an  independent  coil  for  the  flux 
bias.  They  are  protected  independently  by  the  Cryoperm  magnetic  shield  from  an  external  flux 
noise  such  as  the  geomagnetic  field. 

Microwave  pulses  for  the  drive,  signal,  and  pump  ports  are  generated  by  mixing  the  con¬ 
tinuous  microwaves  with  pulses  which  have  independent  IF  frequencies  generated  by  DACs 
(digital  to  analog  converter)  developed  by  Martinis  group  at  UCSB  (<5).  The  pulses  are  applied 
through  the  input  microwave  semi-rigid  cables,  each  with  attenuators  of  42  dB  in  total,  and 
DC-blocks  for  the  drive  and  pump  ports.  For  the  signal  port,  the  microwave  pulses  are  further 
attenuated  by  20  dB,  and  are  input  to  the  resonator  through  a  circulator  to  separate  the  input  and 
reflected  waves.  The  reflected  waves  are  routed  to  PPLO  via  three  circulators  (9-11  GHz)  and 
are  reflected  there  again,  and  are  propagated  through  a  low-pass  (fc  =  12.4  GHz)  and  band-pass 
filters  (9-11  GHz),  two  isolators  (9-11  GHz),  and  the  circulator  (9-11  GHz)  with  a  50  fl  termina¬ 
tion.  Finally,  the  signals  are  amplified  by  a  cryogenic  HEMT  amplifier  and  a  room-temperature 
amplifier  with  a  total  gain  of  ~  66  dB,  and  mixed  with  a  local  oscillator  at  an  I/Q  mixer  down  to 
the  IF  frequency.  The  I  component  of  the  reflected  signals  are  sampled  at  1  GHz/s  by  an  ADC 
(analog  to  digital  converter). 

For  the  impedance-matching  measurement  (Fig.  2A),  the  PPLO  is  kept  off.  Namely,  pump 
pulse  is  off  (the  output  from  the  DAC  in  the  pump  port  is  turned  off)  and  u,’rPO  is  far  detuned 
from  cur  so  that  the  PPLO  acts  as  a  perfect  mirror.  In  other  measurements,  the  PPLO  is  kept  on. 


Input-power  calibration 

To  estimate  the  photon  detection  efficiency  precisely,  calibration  of  the  signal  microwave  power 
level  on  the  sample  chip  is  required.  For  the  calibration,  we  measure  the  reflection  coefficient 
as  a  function  of  the  signal  microwave  frequency  u>s  and  the  drive  power  Pd  and  determine  the 
impedance-matching  points  (Fig.  S3).  Here,  we  use  continuous  microwaves  for  both  the  signal 
and  the  qubit  drive,  and  set  the  drive  frequency  at  cca  =  cuge  —  27r  x  46  MHz.  We  observe 
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Figure  S2:  Experimental  setup  diagram. 
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two  dips  representing  the  impedance  matching,  similarly  to  the  inset  of  Fig.  2A.  In  the  limit  of 
weak  signal  power  and  no  intrinsic  loss  of  the  resonator,  these  dips  are  expected  to  appear  at 
the  same  Pd,  where  the  two  radiative  decay  rates  of  the  A  system  are  balanced,  k3 1  =  k32  and 
k4i  =  k42  (7),  where  k^  is  the  radiative  decay  rate  for  the  |i)  \j)  transition  in  the  impedance- 

matched  A  system.  In  the  actual  system,  however,  the  finite  population  in  the  level  |2)  as  well 
as  the  intrinsic  loss  of  the  resonator  weakens  the  elastic  photon  scattering  from  the  A  system, 
and  the  impedance  matching  occurs  when  the  radiative  decay  rates  are  not  balanced,  k3 \  >  k32 
and  k4i  >  K42  (3).  This  yields  a  difference  in  the  drive  power,  Pd iff,  between  the  two  dips. 
Pdiff  is  sensitive  to  the  input  signal  power:  As  we  increase  the  signal  power,  the  level  |2)  is 
more  populated  and  Pdiff  gets  larger.  Note  that  the  small  Pdig  observed  in  the  inset  of  Fig.  2A  is 
attributed  to  the  intrinsic  loss  of  the  resonator,  since  the  pulsed  signal  field  is  sufficiently  weak 
in  this  measurement. 

We  use  Pud  to  calibrate  the  signal  power  level.  We  determine  the  signal  power  level  which 
reproduces  Pdiff  =  6.0  dB  (Fig.  S3B)  by  the  numerical  simulation,  following  Ref.  (7).  In  the 
numerical  simulation,  we  employ  the  following  parameters  which  are  estimated  by  independent 
measurements:  the  qubit  decay  rate  r/27T  =  0.174  ±  0.012  MHz  (during  this  measurement  T) 
shows  T”1  =  919  ±  62  ns)  and  the  ratio  of  the  external  and  total  decay  rates  of  the  resonator 
photon  Kext/ n  =  0.964  ±  0.003  (for  other  parameters,  see  “Device”  section  of  this  supple¬ 
mentary  material).  As  a  result,  the  signal  power  is  estimated  to  be  Ps  =  —145.28  dBm  at 
maximum  (r/27T  =  0.186  MHz  and  Kext/n  =  0.967)  and  Ps  =  —146.02  dBm  at  minimum 


Pd  (dBm) 


Figure  S3:  (A)  Amplitude  of  the  reflection  coefficient  |r|  of  a  continuous  input  signal  as  a 
function  of  its  frequency  ojs  and  the  drive  power  Pd.  Two  dips  corresponding  to  absorptions  of 
the  input  microwave  due  to  the  impedance  matching  are  observed.  (B)  Cross-sections  of  A  at 
ojs/2tt  =  10.227  GHz  (blue  curve)  and  10.262  GHz  (red  curve).  Difference  in  Pd  between  two 
dips  is  Pdiff  =  6.0  dB. 
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(r / 27r  =  0.162  MHz  and  Kext/ n  =  0.961).  Therefore,  Ps  =  —145.65  ±  0.37  dBm.  This  agrees 
well  with  an  independent  estimation  of  Ps  =  —146.0  dBm  by  taking  into  account  the  total 
losses  in  the  input  port. 


Protocol  for  single  photon  detection 

In  the  main  text  (Fig.  1C),  we  show  the  protocol  for  single  photon  detection.  Here,  we  present 
the  detailed  parameters  of  the  pulses  in  the  protocol. 

The  drive  frequency  is  set  at  uid  =  ujge  —  5u,  where  dco  —  2tt  x  49  MHz  (<  2y)  is  the 
detuning  from  the  qubit  energy,  and  is  fixed  through  all  the  experiments  described  in  the  main 
text.  The  drive  pulse  is  synchronized  with  the  signal  pulse  with  a  Gaussian  envelope.  The 
duration  td  of  the  drive  pulse  is  optimized  as  td  =  1.5 ts  +  50  ns  so  that  the  signal  pulse  is 
completely  covered  by  the  drive  pulse  and  is  efficiently  absorbed  by  the  A  system.  In  order  to 
suppress  unwanted  nonadiabatic  qubit  excitations,  the  rising  and  falling  edges  of  the  drive-pulse 
envelope  are  smoothed  by  Gaussian  function  with  a  full  width  at  half  maximum  (FWHM)  of 
2 trise  =  30  ns  in  the  voltage  amplitude. 

The  readout  pulse  (the  frequency  -ujr(\  =  uT  —  2%  =  2nx  10.187  GHz,  the  length  trd  =  60  ns, 
and  the  mean  photon  number  nrd  ~  10)  is  applied  after  the  delay  of  fdeiayi  =  W2+trise  from  the 
center  of  the  drive  and  signal  pulses.  The  reflected  readout  pulse  works  as  a  locking  signal  for 
the  PPLO  output  phase,  and  the  pump  pulse  (the  frequency  wpump  =  2urd,  the  length  fpump  = 
400  ns,  and  the  power  Ppump  ~  —60  dBm)  is  applied  after  fdelay2  =  40  ns.  The  parametric 
oscillation  signal  with  either  0  or  n  phase  is  output  from  the  PPLO  during  the  application  of  the 
pump  pulse,  and  the  data  acquisition  time  of  ~  100  ns  is  required  to  extract  the  phase. 


Dark  count  in  the  detector 

Figure  S4  shows  the  dark  count  probability  in  the  detector,  which  is  the  click  probability  without 
applying  the  signal  pulse  in  the  pulse  sequence  of  Fig.  1C.  The  dark  count  is  mainly  caused 
by  the  nonadiabatic  qubit  excitation  due  to  the  drive  pulse  and  the  imperfect  initialization.  The 
probability  induced  by  the  latter  factor  is  constant  and  is  measured  to  be  0.008  ±0.001,  while  the 
probability  induced  by  the  former  factor  depends  on  the  power  and  the  length  of  the  drive  pulse 
and  remains  finite  even  with  the  Gaussian  envelope.  We  determine  the  dark  count  probability 
before  and  after  each  measurement  of  Figs.  2B,  2E,  and  Fig.  3  and  subtracted  the  averaged  value 
from  the  measurement  result. 
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Figure  S4:  Dark  count  probability  in  the  detector.  The  data  was  taken  ten  times  each  and 
averaged  before  and  after  the  measurement  in  Fig.  2B.  The  dark  count  probability  including 
the  imperfect  initialization  shows  0.014  ±  0.001  at  Pd  =  —75.5  dBm  where  the  single-photon- 
detection  efficiency  hits  the  maximum.  The  error  bars  represent  the  standard  deviation  in  twenty 
identical  measurements. 

Time  constant  of  an  impedance-matched  A  system 

We  denote  the  overall  decay  rate  of  the  resonator  by  k  and  the  radiative  decay  rate  for  the 
|i)  — *  |  j)  transition  in  the  A  system  by  p, .  Figure  S5  shows  k^/k  as  a  function  of  the  drive 
power  Pd,  calculated  based  on  the  experimental  parameters.  In  the  experiment,  we  choose 
Pd  =  —75.5  dBm  where  the  photon  detection  efficiency  r/  reaches  the  maximum.  At  this  point, 
K41/K  =  0.49.  The  time  constant  of  the  impedance-matched  A  system  for  the  voltage  amplitude 
decay,  ta,  is  estimated  to  be  2/tc  ~  20  ns,  where  k  =  R41  +  R42  ~  27r  x  16  MHz.  The  shortest 
signal  pulse  length  is  34  ns  in  Fig.  3B,  which  is  comparable  with  r\. 


Protocol  for  reset 

In  the  main  text  (Fig.  4A),  we  show  the  reset  protocol  for  the  single  photon  detection.  Here,  we 
describe  how  to  optimize  the  parameters  of  pulses  in  the  protocol. 

We  first  apply  a  tt  pulse  with  the  length  of  6  ns  to  directly  excite  the  qubit  from  the  \g,  0) 
to  the  |e,0)  state.  Then,  we  apply  the  drive  and  reset  pulses  to  induce  the  |2)  — y  |3)  — > 
|1)  transition.  To  find  the  operating  point  to  maximize  the  resetting  efficiency,  we  swept  the 
frequency  curst  of  the  reset  pulse  and  the  drive  power  Pdr.  After  fixing  curst  and  Pdr,  we  optimize 
the  drive  pulse  length  tdr,  and  the  mean  photon  number  in  the  reset  pulse  nrst  to  minimize 
P(  |e)).  Finally,  we  measure  P(|e))  as  a  function  of  curst  and  Pdr  using  the  reset  protocol  with 
optimized  parameters.  Parameters  for  the  readout  and  pump  pulses  are  the  same  as  the  ones  in 
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Figure  S5:  Radiative  decay  rates  of  the  impedance-matched  A  system,  which  are  calculated 
based  on  the  experimental  parameters,  as  a  function  of  the  drive  power.  The  two  relevant  decay 
rates,  /c4i  and  k42  or  Asi  and  k32,  become  identical  at  Pd  =  —75.7  dBm,  where  the  impedance 
matching  takes  place. 


Fig.  1C. 

At  the  optimal  reset  point,  P(|e))  shows  the  minimum  value  of  0.017  ±  0.002,  which  results 
in  twice  larger  occupation  of  the  qubit  excited  state  compared  to  0.008  ±  0.001  obtained  in 
the  initialization  in  the  equilibrium  condition.  This  indicates  the  small  probability  of  unwanted 
nonadiabatic  excitations  by  the  drive  pulse  in  the  reset  protocol. 

It  takes  410  ns  to  reset  the  system  and  208  ns  to  detect  the  single  photon  for  ts  =  85  ns. 
Both  of  the  durations  are  determined  by  the  drive  pulse  widths  including  trise  =  15  ns.  The 
qubit  readout  is  completed  by  accumulating  data  for  100  ns  after  tdeiay2  =  40  ns.  The  period  of 
the  single  photon  detection  including  the  reset  protocol  is  ~  760  ns,  which  allows  the  photon 
counting  rate  of  ~  1.3  MHz. 
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